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ABSTRACT 

Shock jump conditions, i.e., the specification of the downstream parameters of the 
gas in terms of the upstream parameters, are obtained for steady-state, plane shocks 
with oblique magnetic fields and arbitrary flow speeds. This is done by combining the 
continuity of particle number flux and the electromagnetic boundary conditions at the 
shock with the magnetohydrodynamic conservation laws derived from the stress-energy 
tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be 
solved analytically. For mildly relativistic shocks, analytic solutions are obtained for 
isotropic pressure using an approximation for the adiabatic index that is valid in high 
sonic Mach number cases. Examples assuming isotropic pressure illustrate how the 
shock compression ratio depends on the shock speed and obliquity. In the more general 
case of gyrotropic pressure, the jump conditions cannot be solved analytically with- 
out additional assumptions, and the effects of gyrotropic pressure are investigated by 
parameterizing the distribution of pressure parallel and perpendicular to the magnetic 
field. Our numerical solutions reveal that relatively small departures from isotropy 
(e.g., ~ 20%) produce significant changes in the shock compression ratio, r, at all 
shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with 
gyrotropic pressure is obtained. In particular, either dynamically important fields or 
significant pressure anisotropies can incur marked departures from the canonical gas 
dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models 
of particle acceleration in gamma-ray bursts and other environments where relativis- 
tic shocks are inferred. The jump conditions presented apply directly to test-particle 
acceleration, and will facilitate future self-consistent numerical modeling of particle ac- 
celeration at oblique, relativistic shocks; such models include the modification of the 
fluid velocity profile due to the contribution of energetic particles to the momentum 
and energy fluxes. 
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1. INTRODUCTION 

Collisionless shocks are pervasive throughout space and are regularly associated with objects as 
diverse as stellar winds, supernova remnants, galactic and extra-galactic radio jets, and accretion 
onto compact objects. Relativistic shocks, where the shock speed is close to the speed of light, 
may be generated by the most energetic events; for example, pulsar winds, blastwaves in quasars 
and active galactic nuclei (e.g., Blandford k McKee 1977), and in gamma-ray bursts (e.g., Piran 
1999). They may naturally emerge as the evolved products of Poynting flux-driven or matter- 
dominated outflows in the vicinity of compact objects such as neutron stars and black holes. As 
relativistic shocks propagate through space, magnetic fields upstream from the shock, at even small 
angles with respect to the shock normal, are strongly modified by the Lorentz transformation to the 
downstream frame. The downstream magnetic fields are both increased and tilted toward the plane 
of the shock and can have large angles with respect to the shock normal. Hence, most relativistic 
shocks can be expected to see highly oblique magnetic fields and oblique magnetohydrodynamic 
(MHD) jump conditions are required to describe them. 

Relativistic shock jump conditions have been presented in a variety of ways over the years. The 
standard technique for deriving the equations is to set the divergence of the stress-energy tensor 
equal to zero on a thin volume enclosing the shock plane and use Gauss’s theorem to generate the 
jump conditions across the shock. For example, Taub (1948) developed the relativistic form of the 
Rankine-Hugoniot relations, using the stress-energy tensor with velocity expressed in terms of the 
Maxwell-Boltzmann distribution function for a simple gas. de Hoffmann k Teller (1950) presented 
a relativistic MHD treatment of shocks in various orientations and a treatment of oblique shocks 
for the nonrelativistic case, eliminating the electric field by transforming to a frame where the flow 
velocity is parallel to the magnetic field vector (now called the de Hoffmann-Teller frame) . Peacock 
(1981), following Landau k Lifshitz (1959), presented jump conditions without electromagnetic 
fields, and Blandford k McKee (1976), also using the approach of Landau k Lifshitz (1959) and 
Taub (1948), developed a concise set of jump conditions for a simple gas using scalar pressure. 
Webb, Zank, k McKenzie (1987) provided a review of relativistic MHD shocks in ideal, perfectly 
conducting plasmas, and in particular the treatment by Lichnerowicz (1967, 1970), which used this 
approach to develop the relativistic analog of Cabannes’ shock polar (Cabannes 1970), whose origins 
also lie in Landau k Lifshitz (1959). Kirk k Webb (1988) developed hydrodynamic equations using 
a pressure tensor, and Appl k Camenzind (1988) developed relativistic shock equations for MHD 
jets using scalar pressure and magnetic fields with components B z and B# (a parallel field with 
a twist). Ballard k Heavens (1991) derived MHD jump conditions using the stress-energy tensor 
with isotropic pressure and the Maxwell field tensor. By using a Lorentz transformation to the de 
Hoffman-Teller frame, they restricted shock speeds, uq, to uq/c < cos 0 b, where ©b is the angle 
between the shock normal and the upstream magnetic field and c is the speed of light; hence, this 
approach may only be used for mildly relativistic applications. 

All of these approaches assumed that any non-thermal particles encountered by the shock did 
not affect the shock structure, i.e., accelerated particles were treated as test particles, and did 
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not contribute significantly to the MHD flux conservation relations. Moreover, except for Kirk 
& Webb (1988), they confined their analyses to cases of isotropic pressures, a restriction that 
is appropriate to thermal particles or very energetic ones subject to the diffusion approximation 
in the vicinity of non-relativistic shocks. The assumption of pressure isotropy must be relaxed 
when considering the hydrodynamics of relativistic shocks (in particular non-linear ones), since 
their inherent nature imposes anisotropy on the ion and electron distributions: this is due to the 
difficulty particles have streaming against relativistic flows. The computed non-thermal particle 
anisotropies for relativistic shocks are considerable in the shock layer (e.g., Bednarz & Ostrowski 
1998; Kirk et al. 2000), persisting up to arbitrarily high ultrarelativistic energies. Such energetic 
particles eventually relax to isotropy in the fluid frame far downstream, on length scales comparable 
to diffusive ones. However, only a small minority of accelerated particles achieve isotropy in the 
upstream fluid frame of relativistic shocks, due to the rapid convection to the downstream side of the 
flow discontinuity. These isotropized particles are present only when they manage to diffuse more 
than a diffusive mean free path, A, upstream of the shock; their contribution to the flow dynamics 
is therefore dominated by that of the anisotropic particles within a distance A of the shock. For 
non-linear particle acceleration, where the non-thermal ion or electron populations possess a sizable 
fraction of the total energy or momentum fluxes, calculating pressure anisotropy will be critical for 
determining the conservation of these fluxes through the shock transition (e.g., Ellison, Baring, & 
Jones 1996). The jump conditions we develop here can serve as a guide to self-consistent solutions 
of non-linear relativistic shock acceleration problems (Ellison &: Double 2002). 

Here, we extend previous work by deriving a set of fully relativistic MHD jump conditions 
with gyrotropic pressure and oblique magnetic fields. We adopt the gyrotropic case as a specialized 
generalization because it (i) is exactly realized in plane-parallel shocks and is a good approximation 
for oblique shocks where the flow deflection is small, i.e., the field plays a passive role (generally high 
Alfvenic Mach number cases), and (ii) permits a comparatively simple expression of the Rankine- 
Hugoniot jump conditions. The sonic Mach number, Ms, used throughout our paper, refers to the 
shock speed compared to the fast-mode magnetosonic wave speed as described in Kirk & Duffy 
(1999), i.e., Ms = ^3pqUq/(5Po), where p 0 is the upstream mass density and Po is the upstream 
pressure. The Alfven Mach number we refer to here is defined as Ma = -\/47rpon'o/I?o (Ro is the 
upstream magnetic field), regardless of the shock Lorentz factor, 70 = [1 - (uo/c) 2 ] -1 / 2 ; i.e., we use 
these definitions applicable to non-relativistic flows as parameters for the depiction of our results 
at all 70 . 

Our results are not restricted to the de Hoffmann-Teller frame and apply for arbitrary shock 
speeds and arbitrary shock obliquities. We solve these equations and determine the downstream 
state of the gas in terms of the upstream state first for the special case of isotropic pressure, and 
then, by parameterizing the ratio of pressures parallel and perpendicular to the magnetic field, for 
cases of gyrotropic pressure. 

A principal result of this analysis is that either dynamically important magnetic fields or sig- 
nificant pressure anisotropies produce marked departures from the canonical value (e.g., Blandford 
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h McKee 1976; Kirk & Duffy 1999) of r = 3 for the shock compression ratio in an ultrarelativistic 
fluid. The magnetic weakening of ultrarelativistic perpendicular shocks came to prominence in 
the work of Kennel &c Coroniti (1984) on the interaction of the Crab pulsar’s wind with its envi- 
ronment. The similarity of such consequences of fluid anisotropy and low Alfvenic Mach number 
fields, which are also pervasive for trans-relativistic and non-relativistic shocks, has its origin in the 
similar nature of the plasma and electromagnetic contributions to the spatial components of the 
stress-energy tensor. This result may have important implications for the application of first-order 
Fermi shock acceleration theory to gamma-ray bursts and jets in active galaxies. 

In this work, we concentrate on using analytic methods for determining the fluid and elec- 
tromagnetic characteristics of the shock and do not explicitly include first-order Fermi particle 
acceleration. Future work will combine these results with Monte Carlo techniques (e.g., Ellison, 
Baring, & Jones 1996; Ellison & Double 2002) that will allow the modeling of Fermi acceleration 
of particles, including the modification of the shock structure resulting from the backreaction of 
energetic particles on the upstream flow at all pertinent length scales. The jump conditions we 
present here, however, apply directly to test-particle Fermi acceleration in shocks with arbitrary 
speed and obliquity. 


2. DERIVATION OF MHD JUMP CONDITIONS 
2.1. Steady-State, Planar Shock 

Using a Cartesian coordinate system with the -fx-axis pointing towards the downstream direc- 
tion, we consider an infinite, steady-state, plane shock traveling to the left at a speed uo with its 
velocity vector parallel to the normal of the plane of the shock as shown in Figure 1. The upstream 
fluid consists of a thin, nonrelativistic plasma of protons and electrons in thermal equilibrium with 
Tpo — T e o, where T p o (T e o) is the unshocked proton (electron) temperature. A uniform magnetic 
field, Bo, makes an angle ©bo with respect to the x-axis as seen from the upstream plasma frame. 
We keep the field weak enough to insure high Alfven Mach numbers (i.e., M\ > 2.5) and thus to in- 
sure that the magnetic turbulence responsible for scattering the particles is frozen into the plasma. 
The xyz coordinate system is oriented such that there are only two components of magnetic field, 
B x o and B z o, in the upstream frame. The field will remain co-planar in the downstream frame and 
the downstream flow speed will be confined to the x-z plane as well (e.g., Jones & Ellison 1987). 

In the shock frame, the upstream flow is in the +x-direction and is described by the nor- 
malized four- velocity: 0% = 7o(l, /3q , 0, 0). The downstream (i.e., shocked) flow four- velocity is 
@2 — 72(1) /#f , 0,/9f)> where the subscript 0 (2) refers, here and elsewhere, to upstream (down- 
stream) quantities, 0 = u/c, and 7 = (1 - /3 2 ) -1 / 2 is the corresponding Lorentz factor associated 
with the magnitude of the flow three-velocity u. Note that these subscript conventions follow those 
used in Ellison, Berezhko, & Baring (2000), where the subscript 1 is reserved for positions infinites- 
imally upstream of the subshock discontinuity, admitting the possibility of flow and field gradients 
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Fig. 1. — Schematic diagram of a plane shock, viewed from the shock frame, showing oblique 
magnetic fields and a Gaussian volume over which the divergence of the stress-energy tensor is 
integrated. The angles and field strengths are measured in the local plasma frame and, in all of 
our examples, we take the upstream flow to be parallel to the shock normal. The y-axis is directed 
into the page. 
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upstream of the subshock. Here the use of subscript 1 is redundant. In addition, greek upper and 
lower indices refer to spacetime components 0-3, while roman indices refer to space components 
1-3. 

The set of equations connecting the upstream and downstream regions of a shock consist of the 
continuity of particle number flux (for conserved particles), momentum and energy flux conserva- 
tion, plus electromagnetic boundary conditions at the shock interface, and the equation of state. 
The various parameters that define the state of the plasma, such as pressure and magnetic field, 
are determined in the plasma frame and must be Lorentz transformed to the shock frame where 
the jump conditions apply. We assume that the electric field is zero in the local plasma rest frame. 
In general, the six jump conditions plus the equation of state cannot be solved analytically because 
the adiabatic index (i.e., the ratio of specific heats, whose value varies smoothly between the nonrel- 
ativistic and ultrarelativistic limits) is a function of the downstream plasma parameters, creating 
an inherently nonlinear problem (e.g., Ellison & Reynolds 1991). Even with the assumption of 
gyrotropic pressure, there are more unknowns than there are equations; however, with additional 
assumptions, approximate analytic solutions may be obtained. 


2.2. Transformation Properties of the Stress-Energy Tensor 


The stress-energy tensor, T v “', describes the matter and electromagnetic momentum and energy 
content of a medium at any given point in space-time. Continuity across the shock is established by 
setting the appropriate divergence (or covariant derivative T 1 * 1 '- v ) of the stress-energy tensor equal 
to zero, namely for directions locally normal to the shock. Following standard expositions such as 
Tolman (1934) and Weinberg (1972), the total stress-energy tensor, is expressed as the sum 
of fluid and electromagnetic parts, i.e., 

J tot ~ 1 fluid ^ 1 EM ' t 1 / 

The divergence is integrated over a thin volume containing the shock plane as shown in Figure 1. 
Application of Gauss’s theorem then yields the energy and momentum flux conditions across the 
plane of the shock by using T^riy = 0. Accordingly, T QL/ n u yields the conservation of energy flux, 
T lu n u yields the ^-contribution to momentum flux conservation in the ^-direction, T Zu n v yields 
the ^-contribution to momentum flux conservation in the ^-direction, and n„ = (0, 1,0, 0) is the 
unit four-vector along the x-axis in the reference frame of the shock. The Einstein summation 
convention is adopted here and elsewhere in this paper. 


The components of the fluid and electromagnetic tensors are defined in the local plasma frame 
and are subsequently Lorentz-transformed to the shock frame where the flux conservation conditions 
apply, i.e., 

I n at 

( 2 ) 


rnpv — A ^ 


qna/3 


4 . 

+ i EM 


,] . 


where the subscript n (s) refers to the plasma (shock) frame with the x-axis oriented normal to the 
shock in each case. Since the flow speeds in our model may have two space components, in the x- 
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and z-directions, the Lorentz transformation is 


( 


AM = 


iPx 


7 Pz \ 




o o 


iPz -PA 0 1 y 


( 3 ) 


where the /? components and 7 are as previously defined in Section 2 . 1 . The system defined in 
Figure 1 is invariant under translations in the y-direction, so that conservation laws along the 
y-axis are trivially satisfied. 


2.3. The Fluid Tensor and Equation of State 


The fluid tensor will be constrained to the gyrotropic case in the local fluid frame; i.e., pressure 
can have one value parallel to the magnetic field and a different value perpendicular to the magnetic 
field (with symmetry about the magnetic field vector). This gives the diagonal stress-energy tensor 


2 fl u id,B = 


e 0 0 

0 P„ 0 

0 0 Pl 

0 0 0 


0 \ 
0 

0 


P±J 


( 4 ) 


where e is the total energy density, P|| (Pj_) is the pressure parallel (perpendicular) to the magnetic 
field 4 , and the subscript B refers to the magnetic axis in the plasma frame. We obtain the fluid 
tensor in the xyz plasma frame (labeled with subscript n referring to the shock normal), Tfl U id,n, 
with a rotation about the y-axis, 


Tfluid,n — Tfluid.B 1 , 

where 

1 0 0 0 \ 

_ 0 cos Ob 0 sin 0 b 

“0010 

0 -sin©B 0 cosGb / 


( 5 ) 

( 6 ) 


4 Sucli a decomposition into projections was used, for example, by Chew, Goldberger & Low (1956). 
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The resulting tensor in the xyz plasma frame is 


2fluid,n — 


_ 0 P|| cos 2 ©b + P± sin 2 ©b 0 (Pi — F||)sin0BCOS0B 


0 (Pi — P|j) sin 0 b cos 0 b 0 Pi cos 2 0 b + P|| sin 2 ©b 


with the T l i space components corresponding to the 3-dimensional pressure tensor presented by 
Ellison, Baring, & Jones (1996) for nonrelativistic shocks, once a typographical error is corrected. 

Renaming the components of the fluid tensor as 


Pfluid.n — 


0 0 


with Pxz = Pzx, we have the identities 


Pit = Prx — Pc 


Pi — Prx d" Px 


sin©B 

’ COS 0B ’ 

cosQb 
'sin©e ’ 


P zz — P xx + 2 P zz cot (20b) • (If) 

The T 00 component, e, is the total energy density in the rest or plasma frame. The other com- 
ponents, Pij, are defined by Tolman (1934) as the “absolute stress” components in the proper 
frame. Pij is the pressure parallel to the i-axis exerted on a unit area normal to the j- axis. Hence, 
the diagonal components can be considered a pressure, but the off-diagonal components are shear 
stresses. The isotropic scalar pressure, P = 7V(Pfl u id,B)/3, is a Lorentz invariant. The fluid tensor, 
in general, changes its appearance significantly under a Lorentz transformation via the mixing of 
components: for example, it picks up momentum flux components in reference frames moving with 
respect to the proper frame (Tolman 1934). 

Using the above expressions, we can derive an adiabatic equation of state. Starting with the 
spatial portion of the fluid tensor: 


Pfluid.B = 


the total energy density can be written as 


P,l 0 0 

0 Pi 0 
0 0 Pi 


= JV(Pfluid B) 2 

3(r — l) + 9 ’ 
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where Tr(Va u id,B) is the trace of the pressure tensor, T is the adiabatic index, and pc 2 is the rest 
mass energy density. Using equations (9) and (10), 


= 3fF=Ty (p » +2P - l) = r^r [ p “ + % (2cot 0B _ tan0B )] ■ < 14 > 

In terms of the magnetic field, where tan 0 b = B z /B x , the adiabatic, gyrotropic equation of state 
becomes 

1 r D/D D\l 

(15) 


_ 1 p . Pxz ( o B Z \ i nr 2 

- r _i[^xx+ 3 \2 Bz Bx ) j+pc . 


While T is well defined in the nonrelativistic and fully relativistic limits, in mildly relativistic shocks 
T depends on an unknown relation between e and P which we will address in a later section. 


2.4. The Electromagnetic Tensor 


The electric and magnetic field components of the general electromagnetic tensor in the plasma 
frame are given by Tolman (1934) as 


^EM,n — 


where the Q’s are the Maxwell stresses defined as 

E 2 4- B 2 j-,2 d 2 



e 2 + b 2 
2 

(E x B) x 

(E x B) y 

(E x B), ^ 

_1_ 

(E x B) x 

Qxx 

Qxy 

Qxz 

4ir 

(E x B) y 

Qyx 

Qyy 

Qyz 


{ (E x B), 

Qzx 

Qzy 

Qzz J 


Qii — 


— E± — B+ and Qij == — ( E{Ej + B{Bj) . 


(16) 


(17) 


Note that here the suffix n denotes an axis orientation along the shock normal, as it does for the 
fluid tensor. Since the electric fields in the plasma frame are negligible, this simplifies to 


^EM,n — 4^ 


/ B^_ 

2 

0 
0 
0 


0 

2 

0 

-B Z B X 


0 

—B X B Z 


0 
0 

0 


\ 


E - £ 


(18) 


Observe that this electromagnetic contribution to the stress-energy tensor resembles the structure 
of that for the gyrotropic fluid in equation (8), i.e., the presence of laminar fields should mimic 
anisotropic pressures in terms of their effect on the flow dynamics. A noticeable difference, however, 
is that while the magnetic field can exhibit tension and can therefore generate negative diagonal 
components for T #u ', depending on the orientation of the field, the corresponding diagonal pressure 
components are always positive definite. 
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2.5. Flux Conservation Relations 

As discussed above, the energy and momentum conservation relations in the shock frame can 
be derived by applying T^ v n v = 0 to equation (2), individually on the Lorentz-transformed fluid 
and electromagnetic tensors. The conservation of energy flux derives from 

Z& s n v = Ttd.s n„ + s n„ = 0 . (19) 

The fluid contribution to energy flux conservation is 

f£" = = 7 2 A(« + P..) - 7(7 - - P ») + 

7 [(2 7 -lW|+^]|fP,. , (20) 

while the electromagnetic contribution is 

if n M = T^ )S n„ = [( 7 - 1 )MlBl + ( 7 $ + PDPzB* - {(2 7 - l)/3 2 + /3 2 2 }(3 Z B X B Z ] . (21) 


The conservation of momentum flux derives from 

Ttot,s n u = ^flu id)S n„ + Tem,s n„ = 0 . 


( 22 ) 


The ^-component of the transformed fluid tensor contributing to the conservation of momentum 
flux is 


F™ d = T^ d)S n v = 7 2 ^(e + P xx ) + P xx - (7 - l) 2 ^(Pxx - P«) + 


2^ 2 


2( 7 -l)(7/3 2 + /? 2 )£#Px Z , 


0* 


and the z-component is 

pfluid qn3iy _ ^,2 


2 \§xPz 


F*r = 2fluid, s », = 7 2 &&(e + ^x) - (7 - 1)(« + 7 «)“P(JV. - P«) + 


7 + 2 ( 7 - 1 ) 2^1 


(23) 


(24) 


The x-component of the electromagnetic contribution to the conservation of momentum flux is 


77EM rp\v „ 

2px — EM,s n v 


= j^B 2 + ^ [( 7 Pl+Pl? - (7 - i)VM] (B. 2 - Bl) 


2?r 


(7 - 1)(7, 


(25) 
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and the z-component is 


F™ = ^EMfi n f — sjAAS 2 + Sfh - 1) 2 (/*J - 


2n 


(7 


2 

z 

d2 o2 


bD- 


1 ) 2 ^B x B z - £ Bx B z 


(26) 


Clearly, parallels between the various components of the fluid and electromagnetic stress tensors 
can be drawn by inspection of these forms for the fluxes. 


In all cases where the Alfven Mach number is greater than a few, the downstream flow velocity 
deviates only slightly from the shock normal direction so (3 Z ^ fix- This allows a first-order 
approximation in /3 Z and the above equations become: 


ifn id « 7 2 Ar( e + Pxx) + 7(27 — l)PzPxz , 


riEM 
r en 



7(27 - 1) 

47T 


PzB x B z 




(27) 


for the energy flux contributions, and 


riflllid 

• px 


riEM 

~px 




lPx( e + Pxx) + Pxx + 27(7 - 1 )fPxz , 

p2 , 7_( d2 — B 2 ) — — 11 §£ R R 

8n Px±l + 8n [LS * ^ 2 n P x xtSz ’ 


(28) 


and 

« 1 ! UA( t + P„) + -iP„ - h - - P „) , 

for the x— and z — components of momentum flux, respectively. As we show below, the approxi- 
mation, /?2 /?*, becomes progressively better as the shock Lorentz factor increases but, in fact, 

equations (27-29) provide an excellent approximation at all Lorentz factors for virtually the entire 
parameter regime we have considered in this paper. Unless Mach numbers less than a few and/or 
extreme anisotropies are considered, equations (27-29) yield solutions within one part in 10 4 to 
those obtained with equations (20-26) and are much easier to solve. 



2.6. Jump Conditions 

The jump conditions consist of the energy and momentum flux conservation relations, the 
particle flux continuity, and the boundary conditions on the magnetic field. The conservation of 
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particle number flux 5 is 



where the brackets provide an abbreviation for 


72^2/5*2 - 7ono/?x0 = 0 . 


(30) 


(31) 


This jump condition, as well as the ones that follow, are written in the shock frame and, as always, 
the subscript 0 (2) refers to upstream (downstream) quantities. The remaining jump conditions 
are: 


[ rpfiuid i dEM 
r en ' r en 

[ rifluid , rEM 
r px ' r px 


rifluid 

f p* 


, rp EM 
+ *pz 


o, 
0, 
0 . 


(32) 


Adding the steady-state conditions on the magnetic field, namely that V • B = 0: 

i 2 

B x \ = 0 , 


(33) 


and also that V x E = 0: 


"YiPzBx P x Bz) 


= 0 , 


(34) 


completes the set of six jump conditions. In the limit of (7 — 1) 1 and for isotropic pressures, 

the above expressions reproduce the standard continuity conditions at non-relativistic shocks (e.g., 
see p. 117 of Boyd & Sanderson 1969). 


At this point there are eight unknown downstream quantities (/3 i2) Pz 2 , B x 2 , B Z 2 , P X x 2 , Pxz2, 
e 2 , and 712) and only six equations. If isotropic pressure is assumed, P xx = P zz = P and P xz — 0, 
leaving six equations and seven unknowns. To obtain a closed set of equations for isotropic pressure, 
an assumed equation of state (e.g., equation 15) is added to the analysis. Successive elimination 
of variables then generally leads to a 7th-order equation in the compression ratio r = /3 x0 / /3 i2 
with lengthy algebraic expressions for its coefficients (e.g., Webb, Zank, & McKenzie 1987; Appl 
& Camenzind 1988). The forms of the coefficients depend on the assumed equation of state and 
do not simplify easily. Here we perform some of the simpler algebraic eliminations and then use a 
Newton-Raphson technique to iteratively solve two simultaneous equations. 


In Section 3.1 below, we derive an approximate expression for the downstream adiabatic index, 
r 2 , for cold upstream plasmas. For anisotropic cases, further microphysical information is required, 


5 We assume there is no pair creation nor annihilation. 
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which is generally only accessible using computer simulations; in the gyrotropic approximation, we 
parameterize pressure anisotropy in Section 3.2 to provide insight into global characteristics of the 
jump conditions. 


3. RESULTS 


3.1. Isotropic Pressure, High Sonic Mach Number Cases 

Oblique shock jump conditions cannot, in general, be solved analytically even for isotropic 
pressure because the downstream adiabatic index, T 2 , depends on the total downstream energy 
density and the components of the pressure tensor (or scalar pressure), which are not known before 
the solution is obtained. The problem is inherently nonlinear except in the nonrelativistic and 
ultra-relativistic limits where I 2 = 5/3 and 4/3, respectively. Furthermore, the gyrotropic pressure 
components are determined by the physics of the model and do not easily lend themselves to analytic 
interpretation, although Kirk & Webb (1988) provided equations based on a power-law distribution 
in momentum for the pressure tensor components in the special case of a parallel relativistic shock 
with test particle first-order Fermi shock acceleration. 

An excellent approximation can be obtained in the absence of efficient particle acceleration if 
u t h -C uo, where v t h is the thermal speed of the unshocked plasma. In this case, which corresponds 
to high sonic Mach numbers, upon scattering in the downstream frame all particles have 

7rel =* (1 - ^rel)“ 1/2 > ( 35 ) 


where 

O _ ^rel _ PO -P2 
Pre] c 

is the relative /3 between the converging plasma frames. Prom kinetic theory, the isotropic pressure 
is 



P = 


n 

3 <P-V>, 


(37) 


where n is the particle number density, and p and v are the particle momentum and velocity, 
respectively. Then, with our approximation for particle velocity and using the isotropic version of 
equation (13), 


r 2 


, -I (l/3)Prel^rel , . 7rel/?rel , i 

- PC 2 (7rel ~ 1 )mc 2 3(7 re l - 1) 


(38) 


or, 


■n _ 4 7rel + 1 


(39) 


This essentially kinematic approximation, which is operable only if diffusive transport of particles 
from downstream to upstream contributes insignificantly to the momentum and energy fluxes, per- 
mits a direct numerical solution for isotropic pressure, arbitrary obliquity (as long as the upstream 
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Fig. 2.— Compression ratio, r, 0B2, and ©u 2 versus 0 O 7o- The three panels on the left have 
Ma > 100 and the sonic Mach numbers as shown. The three right-hand panels have Ms > 100 
with the Alfven Mach numbers as shown. In all cases, ©bo = 5°. Note that in the top right-hand 
panel, the Ma = 6 and 60 curves are nearly identical. 
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Alfvenic Mach number, Ma, is high), and arbitrary flow speed (see, for example, Kirk 1988; Gallant 
2002, for alternative forms for I^). Note that equation (39) provides an upper limit to the adiabatic 
index because any particles accelerated by the shock would tend to raise the average Lorentz factor 
and cause the adiabatic index to decrease. 

In Figures 2 and 3 we show results for the compression ratio, r, and the downstream angles, 
032 and @ u 2 , as a function of /3o7o, for two extreme upstream magnetic field angles, 0bo, and 
various sonic and Alfven Mach numbers. The magnitude of the downstream field is 

£2 = [^0 + 1oV-1)B*J 1/2 - (40) 

There are a number of important characteristics of these results. First, the jump conditions 
map smoothly from fully nonrelativistic to ultrarelativistic shock speeds and obtain the canonical 
values for the compression ratio r = P x o/Px2 = 4 for high Mach number, nonrelativistic shocks, and 
r = 3 for high Mach number ultrarelativistic shocks. For /?o 7 o < 1 and regardless of the obliquity, 
a low Ms results in a weaker shock with smaller r, as expected (similar behavior is exhibited in 
Figure 1 of Appl k Camenzind 1988). For /?o 7 o £ 10, the sonic Mach number Ms has little influence 
on the results until it becomes very low, i.e. 70 M| ~ 1 . Since our definition of Ms is inherently 
non-relativistic, perceptible changes to the fluid dynamics arise only when the pressure Pq becomes 
comparable to the relativistic ram pressure joPqPoc 2 ] this domain is exhibited in the parallel fluid 
shock jump conditions explored by Taub (1948). 

Figure 2 with 0bo = 5° and Figure 3 with ©bo = 85° illustrate an important characteristic of 
relativistic shocks. For all upstream field obliquities other than ©bo = 0, the downstream magnetic 
field angle shifts towards ©32 = 90° as the shock Lorentz factor increases, indicating the importance 
of addressing oblique fields when treating acceleration at highly relativistic shocks. The transition 
criterion is directly obtainable from equation (34) (since © U 2 is, in general, very small). One 
quickly arrives at /?o7o tan ©bo ^ 1 being the necessary condition to render tan ©32 > 1, so that 
©B 2 -> 90°. Furthermore, in the left hand panels of Figures 2 and 3, the angle the downstream 
flow makes with the shock normal, © U 2 , is small at all 70 A) (note the logarithmic scale for © U 2 ), 
consistent with the assumption that /3 Z C /3 X . This is a consequence of the passive magnetic field 
corresponding to Ma 1. Very different behavior is exhibited when the field becomes dynamically 
important, as is evident in the upper right-hand panels of Figures 2 and 3: a low Ma (i.e., high 
Bo) can drastically lower r, even at ultrarelativistic speeds (see Kennel k Coroniti 1984). 

In Figure 4 we show T 2 for the two extreme sonic Mach number cases from Figure 2. Our 
approximation for 1 ^ depends on uq and r, but turns out to be quite insensitive to r, at least in 
the range above r ~ 2.7. The generated by equation (39) is essentially the same as the low 
temperature solutions presented in Figure 2b of Heavens k Drury (1988) for an e + e“ plasma. As 
noted above, if Fermi acceleration is permitted to occur, I^ will approach 4/3 at lower /?o 7 o due to 
the contribution of energetic particles. 


The variation of r with shock speed we show here closely resembles the low sonic Mach number 
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solutions depicted in Figure 4 of Fujimura & Kennel (1979), who considered jump conditions in 
parallel (i.e., — >■ oo) trans-relativistic shocks using an isotropic Jiittner-Synge equation of state 

(Synge 1957). The principal effect of upstream heating is to lower the compression ratio and 
weaken the shock when the ratio of the particle pressure to the rest mass energy exceeds the square 
of the shock four-velocity. The compression ratio clearly becomes insensitive to the plasma heating 
in ultrarelativistic, parallel shocks since To is always 4/3 and the solution is /3o/?2 = 1/3 (e.g., 
Blandford & McKee 1976). An array of possible jump conditions is admitted when an extension 
to multi-component plasmas is explored, such as in Peacock (1981), Kirk (1988), Heavens & Drury 
(1988) and Ballard & Heavens (1991), where the thermal interplay of ions and electrons on the shock 
dynamics in the trans-relativistic regime can be encapsulated using the two adiabatic shock index 
parameters To and T 2 - Such a parameterization implies that a variety of equations of state can be 
accommodated within the formalism presented here. Notwithstanding, when particle acceleration 
is considered, thermal equations of state become inappropriate and simulation results become more 
important (Ellison & Reynolds 1991). 
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Fig. 5. — Compression ratio versus Alfven Mach number for various upstream magnetic obliquities, 
©bo- The values of r are obtained in the ultrarelativistic limit with a = 1, 
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Figure 5 shows how r varies as a function of Ma in the ultrarelativistic limit for various values of 
0bo and with isotropic pressure. The compression ratio depends strongly on the upstream obliquity 
and can drop well below 3 for low M\. The lower right-hand panels of Figures 2 and 3 also show 
that a low Ma produces magnetic stresses that cause the downstream flow to deflect from the 
shock normal direction, causing 0 U 2 to vary inversely as M\. Despite the fact that 0 U 2 ~ 20° for 
Ms = 100 and Ma = 2.6, the approximate equations (27)-(29) give essentially identical results as 
the complete equations (20-26). The impact of either Mach number on 0b2 is relatively small. 

The most important consequence of a large magnetic field energy density is that it lowers the 
compression ratio in oblique shocks (the field is dynamically passive in parallel ones with 0bo = 0), 
even when 70 > 1. The effect is contained in Eq. (4.11) of Kennel k Coroniti (1984), which specifies 
the jump condition or downstream flow four-velocity for an ultrarelativistic perpendicular MHD 
shock. Algebraically simplifying their formula, and expressing it in terms of the three-velocity 
compression ratio r and Alfvenic Mach number Ma used here, leads to the form (for (3 x o = /?o ~ 1) 

r = « ^{v' M A + l 6M A + 16-2-Ml} . (41) 

where the ratio a of the magnetic plus electric energy flux to the particle energy flux that is used 
by Kennel k Coroniti (1984) is given by a w 1/M| = £$/(47rpoUo)- ^his result applies to 70 > 1 
regimes, and is reproduced in the numerical results depicted in the top right hand panel of Fig 3 
and also the 89° curve of Fig 5. 

The weakening of nonrelativistic shocks at low Alfven Mach numbers is widely understood. 
Such an effect is suggested for shocks with /3 x o £ 0.8 in Ballard k Heavens (1991) and is also 
somewhat apparent for /3 X 0 < 0.97 in Appl k Camenzind (1988). Kirk k Duffy (1999) show r 
as a function of Ma through the trans-relativistic regime for ©bo = 45°. The mildly relativistic 
regime was appropriate for shocked jets in active galactic nuclei, the main application of relativistic 
MHD shock analyses over a decade ago. Here, it is evident that this lowering of r by dynamically 
important magnetic fields persists up to arbitrarily high 70, a result obtained by Kennel k Coroniti 
(1984), who applied relativistic MHD to the consideration of perpendicular pulsar termination 
shocks. The more recent association of ultrarelativistic shocks with cosmic gamma-ray bursts 
further motivates the extension to the 70 > 100 regime. Accordingly, this magnetic weakening 
of the shock has profound implications for the interpretation of gamma-ray burst spectra and 
associated emission mechanisms. 

The origin of this effect can be attributed to the anisotropic and intrinsically relativistic nature 
of the field structure. In the ultrarelativistic regime, the equation of state of a quasi-isotropic, 
turbulent field structure replicates that of the familiar ultrarelativistic gas with T = 4/3. However, 
if the field is laminar and oblique, the stress-energy tensor in equation (18) exhibits anisotropic 
stresses. This anisotropy influences the flux conservation relations for dynamically important fields 
when B z o ^ 0; moreover it should mimic to some extent the behavior anticipated from anisotropies 
due to the gas contribution. Such a parallel obviously emerges from results presented in the next 
Section. 
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Fig. 6. — The top three panels show the compression ratio, r, versus /?o7o for shocks with varying 
downstream anisotropy, a = Pj_/P \ |, and obliquity, ©bo- In all panels, the solid curves have a = 1, 
the dashed curves have a = 0.8, and the dotted curves have a = 1.2 downstream. The upstream 
plasma is taken to be isotropic (i.e., a = 1) in all cases, but the upstream a has virtually no 
influence on the solutions for the parameter regime we consider. The anisotropy has little effect on 
©B2> but does cause modifications in 0 U 2 at mildly relativistic shock speeds, as indicated in the 
bottom panel which shows the angle the downstream flow makes with the shock normal for the 
particular case 0bo = 10°. These are high Mach number examples (Ms = > 100); the effect 

on r at lower Mach numbers is shown in Figure 7. 
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3.2. Gyrotropic Pressure 

For gyrotropic pressure, an additional constraint is needed to close the system of equations and 
obtain an analytical solution; we impose this via 

P± = aP\\ , (42) 

where a is an arbitrary parameter and a = 1 corresponds to isotropic pressure. Equation (42) 
allows us to illustrate the effects of anisotropic pressure, but is not suggested as a model for 
specific acceleration scenarios. In relativistic shocks where the accelerated population contributes 
significantly to the total dynamical pressure, the value of a should deviate significantly from unity 
in both the upstream and downstream regions. Using equation (7), equation (42) then yields 

P xx = -P|| (cos 2 0 b + a sin 2 ©b) , (43) 

Pxz — P\\ (« ~ 1) sin ©b cos ©b , (44) 

and using the fluxes given in equations (20)-(26) or (27)— (29), we have a closed set of equations 
for the jump conditions for shocks with gyrotropic pressure, arbitrary obliquity, and arbitrary flow 
speed. Results for various ©bo’s and a’s are shown in Figure 6 (these results all have Ms = M\ = 
100 but they remain unchanged for larger Mach numbers). The solid curves have a = 1, the dashed 
curves have a — 0.8, and the dotted curves have a = 1.2. In all cases, we have taken the pressure 
in the unshocked gas to be isotropic and a ^ 1 is only applied downstream. 

While our solutions can allow for anisotropic upstream pressure, an upstream a ^ 1 generally 
produces insignificant changes to our results for the parameter regime discussed here. Notable 
exceptions arise when the sonic Mach number is very low, i.e. 70 M| ~ 1. We observe that angular 
distributions at relativistic shocks generated by Monte Carlo simulations (e.g., Bednarz & Ostrowski 
1998; Ellison & Double 2002) and semi-analytic convection-diffusion equation solutions (Kirk et al. 
2000) indicate that for plane-parallel scenarios with ©bo = 0 the accelerated particles are dominated 
by parallel pressure upstream (a < 1) and perpendicular pressure downstream ( a > 1) when near 
the shock. This would suggest a possible weakening of the shock if the non-thermal population 
were to contribute significantly to the dynamics, thereby rendering such a contribution less likely. 
The picture may be much different for oblique and perpendicular shocks. 

The effects of anisotropic pressure on the compression ratio depend strongly on 70 and ©bo- 
When ©bo is small, the downstream angle ©32 is also small at nonrelativistic and mildly relativistic 
shock speeds (top panel of Figure 6). Therefore, at these speeds P xz ~ 0 and P xx ~ Py for any 
a. Since Py = P±/a, the fraction of downstream pressure in P xx is inversely proportional to a 
and since P xx largely determines r, the compression ratio is less than the isotropic value, rj so , for 
a = 0.8 and greater than ri S0 for a = 1.2, as shown in the figure. In contrast, as 70 — > 00, ©32 
approaches 90° for any ©bo > 0 (see Figure 2) and P xx ~ P±_ with P xz again approximately equal 
to zero. Now, the fraction of pressure in P xx will be proportional to a and r < rj S0 for a > 1 
and r > rj S0 for a < 1. The transition where r crosses rj S0 occurs at slower shock speeds as ©bo 



- 22 - 



0 1 2 3 4 5 


Anisotropy, a 

Fig. 7. — Compression ratio versus anisotropy for ultrarelativistic shocks. The top panel shows the 
effects of ©BO) while the bottom panel shows the effects of Mach number keeping ©bo = 89° . Note 
that, depending on« = Pj_/P \ |, the compression ratio r can be either larger or smaller than 3. 
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increases (see the ©bo = 10° panel of Figure 6) until ©bo is large enough (©bo = 60° panel) so no 
transition occurs. While the examples in Figure 6 all have Ms = Ma. = 100, the transition /?o7o 
is independent of Mach number. The dependence of ©b 2 on a is relatively small for the examples 
shown in Figure 6 but © U 2 can change significantly with anisotropy, as shown in the bottom panel 
for the ©B2 = 20° example. 

In Figure 7 we show r as a function of a for various Mach numbers and ©bo’s, all in the 
ultrarelativistic limit. The top panel shows that r is relatively insensitive to the upstream magnetic 
field angle, with the lower panel showing a somewhat greater sensitivity to Mach number. The most 
important aspect of these plots is the fact that r can be either higher or lower than the canonical 
value of 3 if anisotropic pressure is important. This contrasts with the effects of a magnetic field 
(with isotropic pressure), where low Ma’s in the ultrarelativistic limit gave r < 3 only (Figures 3 
and 5). 

3.2.1. Analytic solution for 70 1 

Using our flux conservation relations in the limit 70 > 1, and assuming that the upstream 
pressure is isotropic and the downstream adiabatic index = 4/3, we obtain an analytic solution 
for r valid for any oblique angle and any Ms or M\ greater than one: 

(a -t- l)6or 3 + 2a(wo + bo)r 2 - [2(3a + l)wo + (7a + 3)&oj r + 2(2a + l)(wo + bo) = 0 , (45) 

where wq = eo + Po is the enthalpy density and 60 = B\ z f 4ir. Only the z-component of B appears 
since B x 2 = B x q and this component drops out. Removing the r = 1 root, the solution to the 
remaining quadratic is 

r = 

where q = wq / bo. Solutions in terms of Ma and ©bo can be found using the relation bo = 
Poc 2 sin 2 ©bo/M|. For isotropic pressure, i.e., a = 1, 

r= [(2 + 1 ) 2 + % + 1) 1/2 -(f + 1 ) ’ ( 4? ) 

Equations (47) and (46) reproduce the results obtained with the exact equations shown in Figures 5 
and 7 (bottom panel, solid curve) to a high degree of accuracy. In the limit of Ms > 1 and 
©bo 90°, where the identification q — > M\ « 1/cr can be made, equation (47) reduces to 
equation (41), i.e. the result from Kennel & Coroniti (1984). Once r is determined, the downstream 
z-component of B can be found from equation (40). Observe that equation (46) demonstrably 
exhibits the property dr /da < 0 for all values of q, borne out in the curves of Figure 7. This ability 
of parallel pressure (in this case downstream of the shock) to strengthen the shock contrasts the 
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role of magnetic field tension in weakening a shock, as evidenced in Figure 5 and the analysis of 
Kennel & Coroniti (1984). Note that the analytic expressions that Kirk et al. (2000) derived for 
the angular distribution of accelerated particles at parallel, ultrarelativistic shocks indicates values 
of a > 1 when pitch angle diffusion is operating. 


4. CONCLUSIONS 

We have derived shock jump conditions for arbitrary shock speeds and obliquities. When 
combined with a simple approximation for the ratio of specific heats (equation 39), these equations 
specify the downstream conditions in terms of upstream parameters for isotropic pressure. For the 
case of gyrotropic pressure, an additional arbitrary parameter (equation 42) is required to close the 
set of equations and we have presented a number of solutions where the downstream pressure is 
gyrotropic. The exact equations for energy and momentum conservation are fairly complicated, but 
we have presented simpler, approximate results (equations 27-29) which are extremely accurate in a 
wide parameter regime (i.e., Ms > 5; M\ > 5; 0.8 < a < 1.5) for any shock speed or obliquity. To 
our knowledge, this is the first presentation of oblique shock solutions with gyrotropic pressure that 
continuously span the domains of nonrelativistic, trans-relativistic, and ultrarelativistic shocks of 
arbitrary obliquity. In addition, we have presented an analytic solution for the shock compression 
ratio, r, in the limit of ultrarelativistic shock speeds valid for oblique shocks of any Ms or Ma with 
gyrotropic pressure (equation 46). 

The results presented here assume that no first-order Fermi acceleration occurs, but they ap- 
ply directly to test-particle acceleration where the energy density in accelerated particles is small. 
They also constitute an important ingredient in more complex models of nonlinear particle accel- 
eration. In the test-particle case, the compression ratio r is altered by large magnetic fields and/or 
anisotropic pressures, even at ultrarelativistic speeds. The observation of marked departures from 
the canonical value of r = 3 for the compression ratio of an ultrarelativistic shock, for either dy- 
namically important fields or significant pressure anisotropies, is a major result of this paper. In 
the case of magnetic field influences on the dynamics of relativistic shocks, our results extend the 
conclusions of Kennel & Coroniti (1984) to all shock obliquities. The similarity of such consequences 
of fluid anisotropy and low Alfvenic Mach numbers is a consequence of the similar nature of the 
plasma and electromagnetic contributions to the spatial components of the stress-energy tensor. 

The importance of this result is obvious, since changes in r map directly to changes in the 
power-law index of the accelerated spectrum. This power law is the most important characteristic 
of test-particle Fermi acceleration, and is usually associated with trans-relativistic internal shocks 
in gamma-ray burst (GRB) models (e.g., Rees & Meszaros 1992; Piran 1999). In addition, the outer 
blast wave is an ultrarelativistic shock during most of its active phase, sweeping up and accelerating 
interstellar material. This shock, believed to produce long-lasting afterglows, eventually transitions 
to a non-relativistic phase. Our results are directly applicable to test-particle acceleration models 
of both the internal and external shocks in GRBs, as well as to shocks believed to exist in jets 
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in active galactic nuclei. In future work, we will apply these jump conditions to nonlinear shock 
acceleration models where the accelerated particle population can modify the shock structure. 
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